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Abstract 

We study the singlet-triplet Anderson model (STAM) in which a configuration with a doublet is 
hybridized with another containing a singlet and a triplet, as a minimal model to describe two- level 
quantum dots coupled to two metallic leads in effectively a one-channel fashion. The model has 
a quantum phase transition which separates regions of a doublet and a singlet ground state. The 
limits of integer valence of the STAM (which include a model similar to the underscreened spin-1 
Kondo model) are derived and used to predict the behavior of the conductance through the system 
at both sides of the transition, where it jumps abruptly. At a special quantum critical line, the 
STAM can be mapped to an infinite-?/ ordinary Anderson model (OAM) plus a free spin 1/2. We 
use this mapping to obtain the spectral densities of the STAM as a function of those of the OAM 
at the transition. Using the non-crossing approximation (NCA), we calculate the spectral densities 
and conductance through the system as a function of temperature and bias voltage, and determine 
the changes that take place at the quantum phase transition. The separation of the spectral 
density into a singlet and a triplet part allows us to shed light on the underlying physics and to 
explain a shoulder observed recently in the zero-bias conductance as a function of temperature in 
transport measurements through a single fullerene molecule [Roch N et al, 2008 Nature 453, 633]. 
The structure with three peaks observed in nonequilibrium transport in these experiments is also 
explained. 

PACS numbers: 73.63.Kv, 72.15.Qm, 75.20.Hr, 73.23.Hk 
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I. INTRODUCTION 



For nearly four decades, the Anderson model for magnetic impurities has been the sub- 
ject of intense study in condensed matter physics. Its extension to the lattice (or even the 
impurity model above the so called coherence temperature) describes among others, inter- 
mediate valence systems [l, 2] and heavy fermions 0, 4]. A bosonic version of it has been 
used to describe semiconductor microcavities with strong light-matter interaction The 
Kondo model is derived through a canonical transformation as an integer valence limit of 
the Anderson model Ml. The Kondo effect is also one of the most relevant subjects in 

n 

many-body theory [4J]. A strong resurgence of interest in these many-body phenomena takes 
place in recent years with experimental results in nanoscale systems. Progress in nanotech- 
nology has made it possible to construct nanodevices in which the Kondo physics is clearly 



which constitute 



displayed, for example in systems with one quantum dot (QD) 
ideal systems with a single magnetic impurity in which several parameters can be tuned. 
Scanning tunneling spectroscopy has made it possible to probe the local density of states 
near a single impurity and Fano antiresonances have been observed for several magnetic im- 
purities on metal surfaces [ll, 



12L [13L [ijj] . These antiresonances observed in the differential 



conductance, are a consequence of a dip in the spectral density of conduction states caused 



by the Kondo effect 



m 



16 



17l |. Furthermore, corrals built on the (111) surface of noble 



metals or Cu have been used to project the spectral features of the Fano-Kondo antireso- 



nance to remote places 
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171 ]. The observed Fano line shapes for one magnetic im p urit y on 
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2l|. 



these surfaces have been reproduced by many-body calculations 
The essential physics involved in these nanoscopic systems is well understood in terms of 
the ordinary Anderson model (OAM). In the following, we denote by OAM the simplest 
version of the model, with infinite on-site Coulomb repulsion U, in which a configuration 
with a doublet is hybridized with a singlet. In particular, for systems with one QD with an 
odd number of electrons, the conductance at zero bias is increased below a characteristic 



Kondo temperature Tk as a consequence of the Kondo effect. 



'his is a usual feature of 



lOj or single molecules 221 ] . 



single-electron transistors built with semiconductor QD's jsl, 9, 

In a QD with an even number of electrons, in many cases, the ground state is a singlet 
with all dot levels either doubly occupied with both spin projections or empty. In this 
case, as a gate voltage of either sign is applied, the system goes to a configuration with an 
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odd number of electrons and a doublet ground state. Therefore, the OAM still describes 
the system at intermediate and even electronic occupation. However, in other cases with 
an even number of electrons, due to the strong ferromagnetic (Hund) coupling 23j, it is 
energetically favorable to promote one electron of the occupied level of highest energy to 
the next unoccupied level building a triplet state. When this triplet is well below the other 



states, the system can be described 
exactly solvable by Bethe ansatz {24], 



jy the underscreened spin-1 Kondo model, which is 



251 ] . As a consequence, there is a partial screening of 



;he spin 1 that explains the zero bias Kondo peak observed experimentally in this situation 



26, 
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In fact, in real QD's one expects a second screening channel to be active 



below a characteristic temperature T* suppressing the conductance for bias voltage V or 



temperature T such that eV, kT < kT* 
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301 ] . However, comparison with experiment 
suggests that T* (which depends exponentially on a small coupling constant 30] ) is very 
small, so that one can assume an effective one-channel model for practical purposes 30, |31|. 
When a gate voltage induces a change in the occupation in such a way that the lowest 
state of the isolated dot changes from a triplet to a doublet (or conversely), the appropriate 
model has the form of a generalized Anderson model which has been used to describe valence 



fluctuation between two magnetic configurations 



solved with Bethe ansatz 



35. 
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341 ] . Its impurity version was also 



36[| . In contrast to the OAM (which has a singlet ground 



state), its ground state is a doublet. 

The physical picture becomes more complex and also more interesting when singlet and 
triplet states of the configuration with even number of electrons in the dot lie close in energy 
and no one of them can be neglected (see Fig. [1]). We call the model that describes the 
fluctuations between these states and an odd-particle doublet, the singlet-triplet Anderson 
model (STAM). Again, rigorous results for this model can be borrowed from previous studies 
of intermediate valence systems. Allub and Aligia proposed the model to describe the low 
energy physics of Tm impurities fluctuation between the 4f 12 and 4f 13 configurations in a 
cubic crystal field 37|. Using the numerical renormalization group (NRG) the authors found 
a singlet or a doublet ground state depending on the parameters. Therefore, the system has 
a quantum phase transition when the wave function is forced to evolve continuously between 
these two competing ground states. 

Quantum phase transitions is another topic of great interest in condensed matter physics 
381 ] . Recently, Roch et al. performed several transport measurements through a Ceo QD 
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FIG. 1: Scheme of the lowest lying levels of the STAM: a doublet (a singlet and a triplet) for the 
configuration with an odd (even) number of particles n. 



with even occupancy inserted in a nanoscale constriction [39(. They were able to tune the 
parameters in such a way that a clear manifestation of the above mentioned quantum phase 
transition was observed. The differential conductance dl/dV as a function of temperature 
and bias voltage has been measured at both sides of the transition 39j. On the singlet side, 
a dip in the conductance at V = is observed in agreement with theoretical expectations on 
models similar to the STAM 
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40, 



41| as well as non-equilibrium measurements performed 



in carbon nanotubes 4jJ. On the other side of the transition, dl/dV as a function of V 



shows a structure with three peaks that has not been quan titatively explained yet. We have 

4211. but the source of the observed 



obtained recently a symmetric three-peak structure 
asymmetry remains to be investigated. As the temperature T is decreased, the zero bias 
conductance G(T) first increases, then it shows a shoulder or a plateau and then increases 
again. The authors state that this behavior is not understood and speculate that the increase 
at the lowest temperatures might be due to the opening of another parallel transport mode 
(3^]. This plateau in the conductance was presented by us in a short paper 42]. 

Recently, a comprehensive study of the physics of a two level quantum dot, using NRG 
has been presented 3JJ. This model contains the STAM as a limiting case, when higher 
energy states can be neglected. The advantage of the STAM is that it has fewer parameters 
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and is the minimal model to describe the quantum phase transition when charge fluctua- 
tions are allowed. An important result of the work of Logan et al. is the derivation of an 
extended Friedel-Luttinger sum rule which relates the occupation of the dot with the zero 
bias conductance at very small temperatures 31|. However, in this work, results at finite 
bias were obtained using approximate expressions and the equilibrium spectral density and 
(as previous works) the results do not provide an interpretation of the above mentioned ex- 
perimental findings of Roch et al. Because of the difficulties in exten ding robust techniques 
to the nonequilibrium case (discussed for example in Refs. 42j and 431]). very few studies 



of this problem for finite bias voltage exist 4jJ, |42 ] . 



In this paper, we present several analytical results which shed light on the behavior of the 
conductance near the quantum phase transition. We also present numerical results obtained 
using the non-crossing approximation (NCA), which provide an interpretation of the recent 
experiments of equilibrium and nonequilibrium conductance in Cqq QD's near the quantum 
phase transition 3^]. 

In Section [IT] we explain the model and is application to multilevel QD's. In Section [LTTl 
we derive the integer valence limits of the model by means of canonical transformations, and 
use known results of the ensuing effective models to predict the behavior of the conductance 
at both sides of the quantum phase transition. The self-consistent system of equations of the 
NCA approximation and the expression that gives the current through a system described 
by the model, are presented in Section IIV1 In Section |V] we describe how the STAM for 
a particular set of quantum critical points, can be mapped into an OAM plus a free spin, 
and derive useful results from this mapping. In particular a formula is derived, which allows 
one to calculate the spectral densities of the STAM in terms of those of the OAM. We also 
show that the NCA equations satisfy exact results derived from this mapping. Section IVII 
contains the numerical results obtained with the NCA, and comparison with experiment 
and previous works. Section IVHI is a summary and discussion. 
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II. THE MODEL 



A. The mixed valence impurity 



As originally derived for Tm impurities in a cubic crystal field [37|], the STAM hybridizes 
the lowest states of the 4f 12 configuration, a I\ singlet and T 4 triplet, with a doublet (r 6 or 
T 7 ) of the 4f 13 configuration. The fact that only two neighboring configurations are allowed 
implies that infinite Coulomb repulsion U is assumed. This assumption is taken in all models 
discussed in this paper. Using the notation \ SM), where S is the spin and M its projection, 
we represent the states of the 4f 12 configuration, as 1 00) for the singlet and |1M), (M = — 1, 
or 1) for the triplet. The doublet is denoted by its spin 1/2 projection \a). 
The Hamiltonian can be written in the form 

H = £ a |00>(00| +E t J2 \IM)(1M\ + E d ^ \a)(a\ + H hand + # mix , (1) 

M a 

where i^band is a band of extended states 

Hhand = £k C kcr C k<Ti (2) 

Act 

and -ff m ix is the hybridization. In the following we assume k independent matrix elements. 
Extensions to a more general case is straightforward within the NCA Q, 45]. We can also 
assume full rotational symmetry and then, the form of H m i x is determined by Clebsch-Gordan 
coefficients Calling e£ = Ylk^ka/^^ one obtains 

tfmix = {[y s |t)(oo|-T4(|t)(io| + V2| i)(i-i|)h 

+ [V s | l)(00|+y t (| I)(10| + v / 2| T)(ll|)]q + H.c.}. (3) 

Performing an electron hole-transformation h-[ = — ct, = cj, H mix takes the equivalent 
form 

#mix = {V s (h\\ |) - h\\ T))(00| + V t [{h\\ I) + h\\ t))(io| 

+ v / 2(/i t T |T)(ll| + M|I)(l-l|)] + H.c.}, (4) 
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y,y 



the Bethe ansatz) and the 



361 ] . Thus, the model has 



which is more transparent: the states \a) may be though as having one particle (4f hole 
in the Tm case) and when another particle comes from the band, a localized two-particle 
singlet 1 00) or a component of the triplet |1M) is formed. In any case, it is clear that the 
above transformation allows to treat with the same Hamiltonian the cases in which the 
configuration with the doublet has either one more particle or one less particle than the 
other one. 

We can assume that V s > changing if necessary, the phase of |00). Similarly we assume 
V t > 0. For E t — > +oo, the model reduces to the OAM. For EL — > +oo, the model describes 
valence fluctuations between two magnetic configurations Q, Q] • In both limits, for constant 
density of conduction states, the model is exactly solvable 
ground state is a singlet (doublet) in the first (second) case 
a quantum phase transition as a function of E s — E t . The position of the transition depends 
on the other parameters of the model, leading to a quantum critical surface that can be 
determined calculating the magnetic susceptibility at T — > using numerical renormalization 
group (NRG) [37|]. However, as shown in Section [V] if V t = V s , the transition takes place 
exactly at E x — E t = 0, independently of the value of Ed. In addition, along this quantum 
critical line [46|, the model can be mapped into an OAM plus a free spin 1/2. 

B. The multilevel dot 

In a dot hybridized with two leads, the triplet is formed from the singlet with lowest 
energy by promoting an electron from the highest occupied level, which we denote as a, to 
the highest unoccupied level b. One can restrict to these two levels. The states of these 
two levels are hybridized with the bands of the two leads, left {y = L) and right (v = R) 
described by 

-f^band = e ukhl ka h uka , (5) 

ukcr 

through the following term in the Hamiltonian 

#mix = + VX)hua + H.C.] , (6) 

where h\ a = h [ vka /\fN . We assume V^Vj^ = V^V^, so that only one conduction channel 
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K = C£ K7W/tra 2 + V^) 2 ] 1/2 (v = aor b) (7) 

hybridizes with the dot states. In general, also the orthogonal linear combination of h va 
plays a role and "screens " the remaining doublet ground state when the localized triplet 
is well below the singlet, leading to a singlet ground state [29|, |30J|. However, as mentioned 
in Section [H the characteristic energy scale involved in this second screening T* might be 
exponentially small. As argued before [30], this is likely the case of previous experiments. As 
our results will show (Section IVII) . the one-channel case also describes the recent transport 
measurements in Cqq QD's [39]: the theory in the more general two-channel case 29j, |30j, |47 ] 
predicts that the zero-bias conductance G(T) should decrease at very low temperatures 
and dl/dV should also decrease for the smallest applied bias voltages V in contrast to 
the observations. This indicates that T* is smaller than the smallest temperature in the 
experiments. 

Assuming that the difference between the energies of the levels b and a is larger than 
the hybridization terms, one can retain only the lowest doublet and neglect the singlets that 
contain at least one particle in the state b. Then, performing an electron-hole transformation 
if necessary, the relevant low-energy states of the dot are 



\ a ) = 4|0), |00) = a\a\\0), |11) = b\a\\0) 



|10) 



-L(b\a\ + b\a\)\0),\l-l) 



b\a\\0). 



(8) 



The triplet states should be kept because the ferromagnetic exchange 



lowest of the configuration with even number of particles 23, 



26 



27 



may render them the 



28 



39j (see Fig. 



Restricting the action of the fermion operators to these six states one has 



a\ = |00)(||,6 T T = |ll)(T| + |10)(||/v / 2, 

g{ = -|00)(tUl = |l-l)U| + |10)(t|/V2. (9) 
Replacing Eqs. (C]) and ([9]) in Eq. ([6]) one obtains Eq. (TjJ with 

V s = [(V£f + V^n'\ V t = \{V b L f + (^) 2 ] 1/2 /v / 2. (10) 
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Therefore, at equilibrium, the model takes the same form as Eq. ([I]). The values of E s , E t , 
and Ed are easily determined from the on-site energies and correlations at the QD, including 
Hund exchange [31|] . The nonequilibrium case is discussed in Sections lIV Bl and IIVI 



III. THE INTEGER VALENCE LIMITS 

When | min(E s , E t ) — Ed\ ^> max(V,, Vt), the model is at the integer valence (or "Kondo") 
limit and a Hamiltonian of the exchange type can be derived, using a canonical transforma- 
tion which eliminates H m i x from the Hamiltonian and originates a term quadratic in i^mix 
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|, |48J . Two cases can be distinguished depending on which configuration is favored. 



Odd number of electrons 



For min(E s , E t ) — E^, ^> ma,x(V s ,V t ), the doublet is favored. In this case, the canonical 



transformation is a particular case of that considered in Ref. [2l| for a Cr trimer on Au(lll) 
(a doublet ground state with virtual charge fluctuations to singlets and triplets). Using Eqs. 
(□), P), and (0}, one obtains 



#odd = J2 - £ k h L h ka + (F s - 3F t )hih a + JS ■ S, 
Act 

V 2 

J = A(F S - F t ), F„ = — 3 r, (11) 

V s th n 2{E, n -E d y V ; 

where s = J2 a f3 ha (T af3hp/2 = J2 a /3 c L cr ai3Ci3/'2< is the spin of the conduction electrons at the 
QD and similarly S is the spin of the doublet. This is a Kondo model with potential 
scattering [second term of Eq. ffTT]) ]. It is known from NRG that when J is positive, it 
is a marginally relevant perturbation (it grows with renormalization to J —>■ +oo) leading 



to a singlet ground state [37|, |49[ . Instead for negative J, the exchange term is marginally 
irrelevant (J — > 0) and the ground state is a doublet. Therefore, there is a quantum phase 
transition for J = 0, or F s = F t in terms of the parameters of the original model. 

Moreover, at both sides of the transition, the system is a Fermi liquid, but the phase shift 
at the Fermi energy is 5 = tt/2 for J > 0, but 5 = for J < 0. For the one-channel case 
that we are considering, this means that the zero-bias conductance is maximum (vanishing) 



for J > (J < 0) [29(. Therefore, there is a jump in the conductance at the quantum phase 



transition. This is in agreement with recent NRG results and an analysis of a generalized 



Friedel-Luttinger sum rule 



3l|. 



Note that for V s = V t , the transition is exactly at E s = E t . It has been found by NRG 
that this statement is valid for any occupation of the dot and not only for integer occupation 



371 ] . An analytic argument which supports this numerical result is given in Section IV Bl 



B. Even number of electrons 



For Ed — mm(E s , E t ) 3> max(V^, Vt), if in addition \E S — E t \ 3> maxfV,,, Vt), the highest 
lying .evels between the singlet and the trip.et can be neg.ected, and the resuming enect.ve 
model reduces to one of the cases considered in Ref. [2j. In particular if the triplet is the 

owest in energy one has the underscreened Kondo model with singular Fermi liquid behavior 

25|. 

The effective model H even is reacher when E s lies near E t . In the following we assume 
that both energies lie well below Ed- In order to obtain a more transparent form of H even 
we introduce two fictitious spins 1/2, Si and S 2 , to represent the states of the configuration 
with even number of particles, in terms of the states of these two spins \a1a2) as follows 

|oo) = (| TI) - I IT))/V2, |n) = I TT), |io) = (| t!) + 1 1T»/A |i - i> = I U). (12) 

The resulting effective Hamiltonian is 



2V 2 

#even = V e k c\ a c k(7 + - *— S • (Si + S 2 ) + (E t - E S )(S 1 ■ S 2 ) 

fccr 

( V? V 2 \ 

where n = c\c a and ris = 1 00) (00 1 = 1/4 — Si ■ S 2 . 

For V s = and E t well below E s , H even reduces to the spin 1 underscreened Kondo model 
plus potential scattering. This model has a doublet ground state and singular Fermi liquid 
behavior 25j. A Hamiltonian with the first four terms (the most relevant ones) was studied 



by NRG and found to have a singlet-doublet quantum phase transition which is in general 
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continuous of the Kosterlitz-Thouless type |50j , as the original model [37|. The transition 
is first order only in the particular case in which the fourth term (proportional to Si — S2) 
vanishes 5CJ. This implies either V s — or Vt = in our model. Note that for V s = V t 
and E s = E t , the spin S2 decouples and the spin Si has a usual Kondo interaction with 
the band and is therefore screened. An analysis of the strong coupling fixed point for the 
STAM 37| and for i? eve n without the terms proportional to no 50j indicates that a generic 



feature of the phase with a doublet ground state (and also the quantum critical point) is a 
free spin 1/2 and a Kondo screened spin 1/2 at low temperatures. This implies a phase shift 
5 = 7r/2 and maximum conductance G, as in the phase with a singlet ground state when 
the configuration with odd number of particles is favored (see previous subsection). 

Starting at the transition and decreasing E s , the remaining spin is also screened in a 
second stage and the ground state is a singlet. The full Fermi liquid behavior is restored, 
and one expects 5 = G = 0. Another way to think on this is that the screening of the 
remaining spin leads to a Fano antiresonance with a characteristic energy scale given by the 
second stage Kondo effect (see Section IV Bl) . 

Note that the jump in the conductance at the transition takes place in spite of the fact 
that the transition is continuous (as discussed above) and not first order. The impurity 
contribution to the magnetic susceptibility also jumps at zero temperature at the transition 



23, 



50|. 



These results for the conductance agree with direct calculations using NRG [3l|, |40] 



IV. THE NON-CROSSING APPROXIMATION (NCA) 



A. Representation of the Hamiltonian with slave particles 



The NCA has been used to study the OAM out of equilibrium 52|. To extend this 
formalism to the STAM, we introduce auxiliary bosons, one for the singlet state (s) and 
three for the triplets (t M , M = —1, 0, 1), and auxiliary fermions (f a ) for the doublet, in 
analogy to the SU(N) xSU(M) generalization of the Anderson model 53, |54|. In terms of 
the auxiliary operators, the Hamiltonian takes the form 
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M cr uka 

+ E [(Wt, + vMcuka + H.c] , (14) 



{y — L or R), where Ef = E^ and 



dl = fls, 

4 r = -(ffa + y/2flt-i)/V3, 

d\ { = (/{t + V^/|ti)/V3, (15) 



with the constraint 



^ + E4^ + ^/]/ ff = l. (16) 

The factor l/v 7 ^ in Eqs. (fTS]) was chosen to give a more symmetric form for the NCA 
equations, in particular near the quantum critical line for which the system is exactly solvable 
(see Section IV Dl) . Comparing with Eqs. (JHJ), one can realize that dL = — aj., dj!, = a-f, 
dL- = — dj, = y/2/3b^. The representation of Eq. ( TT4|) was chosen in such a way 
that if the triplet can be neglected (because either E t — > +oo or V[ — Vr — 0), the model 
reduces to the OAM. 

As in Section IIIB| the mixing part of Eq. flT4|) can be put in the form of Eq. ([3]) with 

v s = my + v^n'\ v t = \{vif + (v^/yvz. 



B. Equation for the current 



For the calculation of the current, we consider a multilevel QD with at most six relevant 
states of the Hilbert space, as described in Section III Bl [see Eq. flS])]. For the case of 
proportionate couplings of the relevant two levels {V£Vj^ = V^V^ as we are assuming), Meir 



and Wingreen 



5l| provided an expression for the current in a non-equilibrium situation [Eq. 



(9) of Ref. 51|], which is given by a trace of TG r , where T is given in terms of and G r 
is a matrix of retarded Green's functions. In our case, T and G r are 4x4 matrices in spin 
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and level (a or b) indices. It is easy to see that the product TG r is the same (as it should 
be) in the representation of slave particles used above [Eqs. ffl5|) ]. in which the the index 7] 
refers to s and t instead of a and b (the normalization factors in T and G r cancel). 

r is diagonal in spin index. Within the NCA, the expectation values entering the Green's 



functions decouple into fermion and boson parts |52J and both are diagonal as a consequence 
of SU(2) invariance of the Hamiltonian. This means that G r is diagonal in level index. This 
allows to simplify the resulting expression, which takes the form 

Attp r 

I = — du;J2^P»[h(^) - f R {u>)], (17) 
v 

where p^(co) = — \uiG r dv(T {uj) / ti is the spectral density of & , 

H = Y\ + rl, with K = 2nJ2 IKf% - e k ) (18) 

k 

assumed independent of u within a bandwidth D and zero elsewhere, 

A = 4rK/(Il + T\f < 1 (19) 

(independent of rj) is a parameter that characterizes the asymmetry between left and right 
leads, and f v {u) is the Fermi function with the chemical potential fi u of the corresponding 
lead. 

C. Spectral densities and Green's functions 

The spectral densities of the operators dl a defined by Eqs. f TTol) for given spin, p s d {ui) and 
p d (u)), are determined by convolutions from those of the auxiliary particles p\(uj) with the 
lesser Green's functions Gf(u) (A = s, t, or /) as follows 

p» = I J dJ (G<(u/)p/(u/ + u) + G<{uj' + u)p, n (uj')) , (20) 

Z = J du (Gf(u) + 2Gf(u) + 3G<{u)) , (21) 

where we define Gf{u) as the Fourier transform of (V M (0)tM{t)) (the result is independent of 
M because of SU(2) symmetry) and similarly for Gf(u>) and G"j{uj). The spectral densities 
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of the auxiliary particles are given by the imaginary part of the corresponding retarded 
Green's function as usual: p\(u) = — \mG\{ui) / n . In turn, these Green's functions 



are given in terms of retarded self-energies E^(o>), which as in the case of the OAM 52], 
should be determined selfconsistently. The imaginary parts are given by the following set of 
integral equations 

ImS^H = - J du'T s p f {u')f{u'-u), 
ImS[H = -iy du/T*p f {cS)f(uJ-u,), 

Im^u) = ~\J du' [V s p s {u') + T t p t (cu')]h(u;-u;'), (23) 

where 



f{u) = [rlf L (u)+Tlf R (u)]/T\ 

h(u) = [Tl(l - f L {u)) + Tl(l - f R {u))\/T\ (24) 

and T^, T v (t] = s or t) are given by Eqs. (Tl8l) . The real part of the self energies are obtained 
using Kramers-Kronig relations 

Re£» = -V [ dJ ImS * (a/) . (25) 

Once the retarded self-energies are obtained solving the above system of equations, the 
lesser Green's functions come from the solution of the following integral equations 

G<(u) = \G R (u)\ 2 Z<(u), (26) 



S<(o;) = - J du' T s Gf(u')h(uj' -u), 
E<(u) = — [ dcJT t G5(u')h(u , -u), 

37T J 

= ^f*S [r s C<(o;') + I*G<(a/)]/(" - u>). (27) 
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D. Numerical details 



In the numerical procedure to solve the NCA equations, we have used a set of self- 



adjusting meshes (rather than the fixed one used in Ref. [54]) to describe the spectral 
densities and lesser Green's functions of the auxiliary particles: we have evaluated p\{ui) 
and Gf(u) in the corresponding logarithmic array of discrete frequencies u\, built in each 
iteration in order to have a larger density of points near the corresponding peaks or singu- 
larities of these functions. The procedure guarantees the resolution of the sets of integral 
equations (123]) and (1271) to a high degree of accuracy. To calculate the spectral densities, 
p s d {pj) and p^(w), we have used two different logarithmic meshes centered at the peaks of the 
functions entering Eqs. (T2"U]) . This scheme of numerical resolution allows us to obtain the 
conductance at both sides of the transition, within equilibrium and non equilibrium, and for 
all values of the parameters considered. 

The logarithmic discretization is similar to that used in NRG calculations, where only 
one mesh centered at the Fermi energy fii = and with an arbitrarily large number of 



frequencies near this energy is used 



55] 



V. THE EXACTLY SOLVABLE CASE 



For V s = V t and E s = E t , the model given by Eqs. ([I]), (12]), and (J3j) or (HJ) has additional 
symmetries and is exactly solvable [37] . It has been shown numerically that these equations 



define a quantum critical line (a point for each value of E^/V s ) 46| that separates regions of 



singlet and double ground states [37| . In this Section we provide simple analytical arguments 
to demonstrate these results, map the corresponding spectral densities, and show that the 
NCA is consistent with these results. 



A. The STAM on the quantum critical line 

In analogy to Eq. f]T2l . let us consider a (probably fictitious) system, like the two- level one 
considered in Section [II B] but in which the relevant sing let is |00) = (l/V2){b\a\-b\a\)\0) 
[instead of atat |0) , see Eqs. (jHJ)], and with mixing term 
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#mix = K b thvv + H.c] , (28) 

l'<7 

which does not involve the a\ and a a operators. Clearly, the resulting model [Eq. fl3Tl) with 
E t = E s ] is the OAM with infinite Coulomb repulsion U for the level b (which is exactly 
solvable by Bethe ansatz Q]), while level a within the relevant Hilbert subspace [Eqs. (jHJ) 
with 1 00) replaced as above] reduces to a decoupled spin 1/2: \a) = aj.|0) . Within this 
subspace 

b\ = |ll)( T | + (|10)(I| + |00)a|)/v / 2, 

b\ = |l-l)a| + (|10)(t|-|00)(T|)/v / 2. (29) 
Replacing these equations in Eq. ( |28l) one obtains Eq. (jlj) with 

v s = v t = [{y b L f + {y b R f] l i 2 /V2. (so) 

This shows the equivalence of the STAM for V s = Vt and E s = E t with an OAM plus a free 
doublet. 

B. Effect of singlet-triplet splitting 

Proceeding as above, it is easy to see that for V s = V t = Voam/v^2, but arbitrary E s and 
E t the STAM given by Eqs. ([!]), (|2J), and (jl]), except for an irrelevant constant is mapped 
onto 



H' = -tkh{ a h ka + (E t - E d ) b% + l V OAublK + H.c] 

kcr a a 

+ Ub\b^b\b l + (E t - E s )(S a -S b - 1/4), (31) 

where U —>■ +oo, S a is the spin operator of the spin 1/2 which is free for E t = E s , and 
similarly S 6 = Y. a pb ] a cr a pbp/2. 

For Et — E s , clearly S& is screened as usual in the OAM and the ground state is a doublet 
which is the direct product of a Fermi liquid singlet times the spin state \a). If the state b is 
hybridized with two conducting leads (as above), the physics of the OAM determines that 
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the conductance G = Go sin 2 5, with Go = 2e 2 A/h and 5 = n(^2 a b\b„) j2 [56j. In particular 
when the total number of particles (Y2 a tfjbo) + 1 = 2, one has maximum conductance, in 
agreement with the result discussed in Section IIIIBI 

When the last term of Eq. ( 13TT) is added, one can think S& as representing the spin of 
itinerant electrons in an effective heavy mass Fermi liquid at low energies. In fact, compari- 
son of a mean-field slave-boson treatment with NRG calculations [57| show that this picture 
is qualitatively correct [3] for (^2 a bl r b a ) + 1 ~ 2. Then, as discussed in Section UTT] from 



the physics of the ensuing effective Kondo model [49|, |53] , when E t < E s (ferromagnetic 
coupling) the exchange interaction renormalizes to zero and the ground state continues to 
be a doublet, while for E t > E s a second screening takes place and the ground state is a 
singlet. This results agrees with previous NRG results 33] and confirms that V s = Vj and 
E s = E t corresponds to the quantum critical line. 

If the OAM for E t = E s is in the Kondo regime ((^2 a blb a ) ~ 1), one expects that 
addition of a positive exchange (E t > E,) induces a Fano-Kondo antiresonance, depressing 
the conductance at low temperatures (pM. \^\. while nothing dramatic happens for E t < E s . 
This again agrees with the results of Section IIII1 

While as discussed above, the ground state is a doublet for E t < E s , we remind the reader 
that for realistic two- level systems at low enough temperatures (T < T*, see Section TlIBj) 



a second screening channel should become active 
doublet and a decrease in the conductance 29, 30 



eading to a screening of the remaining 



43]. 



C. Mapping of the spectral densities 



Since the OAM out of equilibrium has been studied before [43j, |52|, |59|, |6Qj, l6l| . results 
for the conductance of the OAM can be extended to the STAM on the quantum critical 
line if one knows how to express the spectral densities p s d {uj) and p^(^) of the operators 
d\ a , d\ a of the STAM [see Eqs f|T5|) ]. which enter the equation for the current ffTTl) . in 
terms of the spectral density Pbijjj) of the operator of the OAM [included in Eq. (|3T|) 
for E t = E s ]. These densities are proportional to the imaginary part of the corresponding 
retarded Green's functions. For example pb{uj) = — ImGKu), where G T b (uj) is the Fourier 
transform of G r b {t) = — iO{t){b a {t)b\ + W a b a {t)) , where 6{t) is the step function. The operators 
d\ a , d\ a of the STAM can be expressed in terms of those of the OAM using Eqs (Tl5l) and 
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the mapping of operators explained in Section IV Al For example 



dti = -^[2fe{a{a T + (b\a\ + &ta{)a;]. (32) 
v 6 

Replacing this into expectation values like (d^(i)dL), one obtains expectation values involv- 
ing six fermion operators. Four of them correspond to a\ and a a operators, which have no 
dynamics (are time independent) and are decoupled from the remaining b\ and b a operators. 
Evaluating the expectation values involving a\. and a a operators, using spin conservation, 
(alofataj.) = 0, and assuming the paramagnetic phase (in particular (a^cia) = 1/2), we 
obtain after some algebra 

(*il(*)<i> = 5<6iM-*)>, (33) 
independently of rj = s or t. Proceeding in the same way for the remaining expectation 
values we finally obtain 

p5(w) = P&(w) = ip6(-«), (34) 

which relates the spectral densities of the STAM to that of the OAM on the quantum critical 
line. 



D. Mapping of the NCA equations 

It is interesting to note that the NCA approach for the STAM described in Section IIVt 
and the corresponding one for the OAM, although they seem to be quite different at first 



glance, can be related and satisfy Eqs. (1341) . The NCA for the OAM [52j makes use of an 
auxiliary boson b and auxiliary fermions f a , describing the electron operator as b a = b*f a . 
In analogy to Eqs. (fTB"]) . the coupling to the right and left leads are described by coupling 
constants and r^. We define the total coupling of the OAM as T = Tr+Tl- The mapping 
between both models described above for V s = V t and E s = E t implies that Y s u = 1/2T V , 
and T l v = 3/2T u . Therefore, T* = 3r s = 3/2r. As a consequence of these relations, the 
equations for the auxiliary bosons s and those for have the same form [see Eqs. (T2"3"|) and 
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(|27|) ]. Moreover the resulting equations for the self-energies take the same form as those of 
the auxiliary fermions f a of the OAM [Eq. (23) of Ref. js^]], but with f u (oo) replaced by 
1 — f v {uS). Similarly, the self-energies of the auxiliary fermions f a of the STAM take the 
same form as those of the auxiliary boson b of the OAM, with the same change as above in 
the Fermi functions f u . This implies the following relations for the auxiliary particles 



p s (u) = pt(u) = Pf(-w), Pf(u) = jog(-w), 

Gf(u>) = Gf(u) = Gj(-u), G<(u) = Gf(-u). (35) 

The density p&(o>) of the real fermion h a of the OAM is given by a convolution [Eq. (21) 
52j] similar to that defining p s d (ui) and p d {ui) [Eq. (BUI) ] with the replacements 



of Ref. 



indicated by Eqs. (|35|) . but with the difference that Z in the denominator is replaced 
by Zqam — j doo (g^{uu) + 2GJ(uj)^. Using Eqs. (121]) and fl35l) . it is easy to see that 
Z = 2Zqam- This lead to Eqs. (13^1) for the relation between spectral densities. We have 
checked it numerically by an independent solution of the NCA equations for both models. 



VI. NUMERICAL RESULTS 



For the numerical solution of the NCA equations, we take V s = Vt, so that independently 
of Ed, the quantum transition occurs at the exactly solvable case E s = E t described above. 
We take T, the coupling of the OAM involved in the mapping described in the previous 
section, as the unit of energy. Therefore V s = 1/2T, and T* = 3/2I\ We take the band 
width of the conduction bands D — 10T. At equilibrium pi = pr and the properties of the 
model depend on E n + pl — E d (77 = s or t) and not separately on the individual parameters. 
Therefore, without loss of generality, we take Ed = Pl = Pr = 0. Out of equilibrium, unless 
otherwise stated, we assume T R = T V L = T r '/2 and Pl + Pr = (as expected for equal 
couplings to left and right leads). 

The choice V s = V t leaves four free parameters: temperature T, bias voltage V which 
determines the difference in chemical potentials Pl~Pr = eV , E t /Y (or E s /V) which controls 
the valence and can be modified by a gate voltage, and finally (E s — E t )/T which controls 
the distance to the quantum critical line j^. This parameter has also been controlled 
experimentally by Roch et al. [39]. A great advantage of taking V s = Vt is that we know 
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exactly where the quantum transition is, while for other ratios V s /Vt, the position of the 



transition has to be determined numerically 



37| . This is very time consuming within the 



NCA because it is required to solve the structure of the spectral densities at low temperatures 
many times near the transition. 

In this paper, we restrict our study to E s ,E t < E d . This means that the configuration 
with even number of particles is favored. This situation corresponds to the most novel 
experimental results, in particular those of Roch et al. for the conductance through Ceo 
QD's near the quantum phase transition [39]. In the following and for the sake of brevity 
we call the "singlet side" of the transition the region of parameters with a singlet ground 
state (E a < E t in our case with V s = V t ), and (to be consistent with Roch et al.) we denote 
by "triplet side" the region E s > E t although as explained above, the spin 1 is partially 
screened and the ground state is a doublet. 



A. The spectral densities 

As shown in Section HVBt the current is proportional to the integral of the following 
weighted average spectral density 

However, as we will show, a study of the singlet p s d (u) and triplet p d (oj) parts of this average 
density contributes significantly to the understanding of the numerical results. As expected 
from Section IV Dl within our numerical accuracy, the three densities, shown in Fi g. J2j 
coincide with the specular image of that of the OAM times a factor 1/2 (Fig. 5 of Ref. [52| ) . 
At equilibrium, they show a peak at the Fermi level. The width of this peak allows to define 
a Kondo temperature Tr-. Under an applied bias voltage, the peak splits in two near the 



corresponding p v as in the OAM 52 ] 



How do the spectral densities evolve as one moves from the quantum critical line? For 
V = 0, this is shown in Fig. [3] on the singlet side of the transition. For small tempera- 
tures, decreasing E s from the quantum critical line (E t — E s ), p s d (uj) displaces to positive 
frequencies, while p d (oS) decreases and displaces its weight to negative frequencies. As a 
consequence, a pseudogap opens in p^ v (cj). A similar pseudogap was found before in studies 
of two-level systems and interpreted as the low temperature part of a two-stage Kondo effect 
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FIG. 2: Spectral densities as a function of frequency for T = 5 x 10 A , Et = E s = —2, and three 
bias voltages 
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FIG. 3: Spectral densities as a function of frequency for V = 0, Et = —2, E s = —2.03 and several 
temperatures. Left: p s d {uf) and p l d {oj), middle: weighted average [see Eq. (f36l) ]. right: p d (u>) at low 
temperatures. 



[40j , along the lines discussed in Sections IIII Bl and IV Bl At very low temperatures (below 
0.005r), a spurious spike appears at the Fermi energy in p%{oj). This is due to a known 
shortcoming of the NCA that takes place when the ground state for zero hybridization is 
non-degenerate, for example under an applied magnetic field 52j. However, as argued in 



Ref. 



52| it is interesting to note that this shortcoming does not affect the calculation of 



thermodynamic properties under a finite applied magnetic field 62] 
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FIG. 4: Spectral densities p d (uj) and p d {oj) as a function of frequency for T = 0.001T, E t 
E s = —2.03 and several bias voltages. 



-2. 



The densities under an applied voltage are shown in Fig. HI The main effect of the bias 
voltage V is to lead to a decrease of the singlet part of the density p s d (uj) and a simultaneous 
increase of the triplet part p d (w) for positive frequencies. In addition, near the average Fermi 
level (/iL+/iij)/2, as V is increased, the pseudogap in the average density of states p^iyj) first 
closes, leading roughly to a single broader peak near (p L + p R )/2, and then, for larger V, this 
oeak in p^i^uj) splits in two near pi and pn as in Fig. [2] and the ordinary Anderson model 
52| | . The right panel of Fig. H] shows p^{oS) with more detail at low energies and including 
another small voltage. This fine structure suggests that as the gate voltage is applied, the 
peaks of p d (uj) split in two, shifted in ±eV/2. Then, naturally, the pseudogap closes when 
eV reaches the difference between the position of the peaks, which is near E t — E s . 
The spurious peak at the Fermi level disappears already for very small bias voltages. 
The spectral densities on the triplet side of the transition (E s > E t ) at equilibrium are 
shown in Fig. [5j In this case, the ground state for V s = V t = is degenerate and no spurious 
peaks appear. Therefore our results are more robust. In contrast to the previous case, p d (w) 
remains peaked at the Fermi energy for low temperatures. This is a consequence of the 
partial Kondo effect, by which the spin 1 at the dot forms a ground state doublet with the 
conduction electrons of both leads, as it is known from the exact solution of the model when 
the singlet can be neglected 
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36] or the spin 1 underscreened Kondo model |24j. |25|. 



The singlet part of the density p s d (uj) displaces to negative frequencies in this case. There- 
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FIG. 5: Spectral densities p d (u>) and p d (io) as a function of frequency for V = 0, E t = —2, 
E s = —1.97 and several temperatures. 

fore a pseudogap also appears in p^(u), but at finite frequencies. Note that at high tem- 
peratures both densities are similar except for a constant factor and the difference in the 
structure of p s d (oj) and p l d {uj) develops at a characteristic temperature of the order of a 
fraction of | E s — Et \ . 

The effect of an applied bias voltage on these densities is shown in Fig. [6j In this case, 
the singlet part of the density increases with voltage at positive frequencies, in contrast to 
the case shown above (Fig. HJ) for the singlet side of the transition. However, the behavior 
of p^ v (o;) near the average Fermi energy is similar as in the above case. The peaks at 
equilibrium first broaden and merge into one for small bias voltage V, and for larger V, this 
peak splits in two centered at energies near pi and /ir. 




B. The equilibrium conductance 



The conductance G(T,V) = dl/dV for V = on the singlet side of the transition is 
shown in Fig. [71 together with the contributions of the singlet [T s p s d (uj) in Eq. fTTTl) ] and 
triple t [r V ^ (q;)| part of the spectral densities. Our result agree with previous ones using 



NRG 



31 



40] and with experiment 



39(. In particular, the increase and decrease of G(T) from 



its maximum value are logarithmic to a good degree of accuracy. At very low temperatures 
(below 0.01 T in Fig. [7J, our result for G(T) increases slightly as the temperature is lowered, 
while one expects a saturation at a value given by the generalized Friedel-Luttinger sum rule 



23 



0.05 




1 1 1 






eV = 0.02 


0.04 




= 0.1 
= 0.2 


0.03 




1 

n 
1 1 


0.02 




1 1 


0.01 




' i i \ 
' 1 V 

■ \ v - 













'}- : i 


-0.5 


0.5 



co/r 



1 1 

r=o.ooir 


- 


i 

! 

I'v 




II 




Pd !! 

1 1 

I 1 




I 1 

^ V: 

'• K 














i 


-0.5 


0.5 


co/r 





0.2 



0.15 



0.1 



-0.05 



FIG. 6: Spectral densities p s d {u) and p l d {uj) as a function of frequency T = O.OOir, E t = —2, 
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FIG. 7: Zero bias conductance G(T, 0) in units of Go = 2e 2 A/h as a function of temperature 
(full line) and contributions from the singlet (dashed dot line) and triplet (dashed line) for Et = 
—3, E s = —3.1. Straight dot lines are guides to the eye. 



3l| (see below). This low-temperature increase is due to the spurious peak that develops in 



p d (ui) as a consequence of the NCA, as explained in Section [VI Al 

According to the generalized Friedel-Luttinger sum rule [311 ]. the conductance G(T, V) 
at zero temperature and without applied bias voltage on the singlet side of the transition is 
given by 

G,(0,0) = G sin 2 (^n odd ) (37) 

where G = 2e 2 A/h and n odd = (E CT W)H) = 1 ~ <|00><00| + J2m |1M)(1M|) is the total 
occupation of the configuration with odd number of particles. For the parameters of Fig. [71 
we obtain n odd = 0.10 ~ and then one expects a low value of G s (00)/G . In fact inserting 
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FIG. 8: Zero bias conductance G(T) as a function of temperature for Et = —3 and several values 
of e = E s — E t (increasing from top to bottom). The circles correspond to Eq. ([38]) . 

this value of n dd in Eq. ( 1371) one obtains G s (0,0)/Go = 0.0245. Our corresponding result 
is near 0.020 (see Fig. [7j). Although it is known that the NCA does not satisfy Fermi liquid 
relationships, the deviation (near 0.005) is not too large for a magnitude that is in general of 
the order of 1. For the particular case of Fig. [7J the deviation is near 10 % of the maximum 
value shown in the figure. 

The evolution of the zero-bias conductance as the system moves from the quantum critical 
line Et = E s to the triplet side of the transition E t < E s is shown in Fig. [HJ At the 
transition, the conductance is the same as that of an OAM 52j, with parameters given 



by the equivalence explained in Sections IV Al and IV CI In particular, we find that the 
conductance is very well described by the empirical curve derived by fitting results of the 
NRG for a spin 1/2 

Ge{T) = r ■ , o1/ . G( °L^ ^ ( 38 ) 



[1 + (2V- -l){T/T K f] s ' 
with s = 0.22. As E t is lowered, removal of degeneracies in the ground state leads to 
a decrease in the conductance at intermediate temperatures. At zero temperature, the 
generalized Friedel-Luttinger sum rule [si|, for the triplet side of the transition gives 



G t (0,0) = G cos 2 (?"odd) • 



(39) 



Since the valence is only slightly increased by a small decrease in E t , one expects that 
G(0, 0) ~ Gq in good agreement with our results. The temperatures reached in our NCA 
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FIG. 9: Zero bias conductance G(T) as a function of temperature (full line) and contributions from 
the singlet (dashed dot dot line) and triplet (dashed line) for E t = —3, E s = —2.9. 

approach for the larger values of e = E s — E t used in Fig. [8] are not low enough to reach 
these high values of the conductance. 

A distinct feature of G(T, 0) on the triplet side of the transition is the developing of a 
bump or a plateau at intermediate temperatures. This is more clearly displayed in Fig. 
[9] where the scale in the conductance has been expanded. This structure has not been 
noticed in previous NRG calculations of the conductance in two-level systems jsj], This 
might be due to insufficient calculations in the appropriate range of temperatures (which 
correspond to rather small values of the conductance). Another possible reason is the loss 



of reso 
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ution of NRG for high-energy features. This shortcoming of NRG is clearly manifest 
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6J, 



65] 



64| in systems of two QD's in which the Kondo resonance is split in two 
The separation of the spectral density p^{uj) which enters the equation for the current ([1 
into singlet and triplet components, as shown in Fig. [9] shows that the bump is due to 
charge excitation involving the singlet component (which are peaked at an energy E s — E t , 
see Fig. [3]) broadened by the temperature. 



This provides an interpretation of the corresponding transport experiments through G 



QD's (Fig. 4 (b) of Ref. 



60 



39j). The comparison suggests that in general, the temperature in 



the experiment could not be lowered significantly after the plateau has been completed, and 
;hat the conductance should continue to increase for decreasing temperatures. Roch et al. 
39l | suggested a different physical picture, fitting the plateau with the empirical Eq. ( |38l) 
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FIG. 10: Differential conductance as a function of bias voltage V for T = 0.01, Et 
values of E s . 



-3 and several 



with a smaller value of Go and speculated that the further increase in G(T, 0) at smaller 
temperatures might be due to the opening another parallel transport mode. The effect of 
a second screening is expected to lead to a decrease of the conductance on general physical 
grounds [29|. Our results also indicate that while Eq. ( l38|) is a good curve fitting for the 
conductance of the OAM, it does not work in the STAM for general values of the parameters. 
An exception is of course, the exactly solvable quantum critical line in which the STAM is 
mapped onto an OAM plus a free spin as described in Section fVT 



C. Conductance as a function of bias voltage 

In Fig. [10] we show the differential conductance G(T, V) = dl/dV as a function of bias 
voltage for small temperatures and for three values of e = E s — E t . One of them (e = 0) 
corresponds to the quantum critical transition and for the other two, the system is either 
on the singlet (e < 0) or triplet (e > 0) side of the transition. The remarkable change of 
behavior at the transition is evident. As expected from the results of Section IVI Al the 
opening of a gap in p™{oS) near the Fermi energy on the singlet side of the transition, leads 
to a dip in G(0, V) at low V. The width of this dip is of the order of |e|. At the bottom 
of the dip, the small value of the conductance G(0, 0) should be given by the generalized 
Friedel-Luttinger sum rule Eq. fl37|) . while the NCA results have a deviation of nearly 20 % 
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FIG. 11: Differential conductance as a function of bias voltage for E% = —3, E s = —2.9 and several 
temperatures. 

of this value, as discussed in Section IVI Al For larger values of — e, the dip becomes wider, 
and the result is similar to the conductance observed in finite chains of an even number of 

Mn atoms on CuN y. 

On the triplet side of the transition, a structure with a central peak at zero bias and 

;hree peaks are more marked at smaller temperatures, 



two lateral maxima is obtained. The 
as we have shown in Fig. 4 of Ref. 42]. This structure also agrees qualitatively with the 
experimental findings in C 60 QD's (Fig. 4 (a) of Ref. [39|). Actually the observed structure 
is asymmetric in contrast to the results shown in Fig. [lOj This is a consequence of our 
assumption of a symmetric voltage drop. In fact the couplings T\ and T V R are usually very 



different for C 60 molecules 



671 ] and this leads to an asymmetric voltage drop. A reasonable 



assumption is that the voltage drop from one lead to the molecule is inversely proportional to 



the coupling to the corresponding lead 60]. In our calculation, the ratio between couplings 
enters through the asymmetry parameter A that controls the magnitude of the current [see 
Section IIVBI and Eqs. (1171) and (119j) ] and the functions ( 1241 that enter the selfconsistency 
equations (see Section HVl) . A calculation for an asymmetric voltage drop is shown in Fig. 
fTTl We have taken \i L = 4eV/5, jiR = —eV/5, and Y r ' R /Y' n L = 4 (assumed independent of 
7] = s or t), keeping the same sums T R + T V L as before. As expected, now the height of 
the lateral peaks is different, with the left peak as the second most intense after the central 
one, in agreement with experiment. The nonmonotonic behavior of G(0, V) can be again 
qualitatively understood from the structure of the spectral densities discussed in Section 
IVI Al Assuming as a first approximation that (P d {oS) does not depend on voltage, it is clear 
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FIG. 12: Full line: differential conductance as a function of bias voltage for E t 
and T = O.OOir. Dashed line: corresponding result taking p^{to) at V = 0. 
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From Eq. ifPTj) that the conductance at zero temperature G(0, V) would proportional to the 
average of p^{^) in a window of uo of width eV around the Fermi energy. Since p^{ui) is 
peaked at the Fermi energy (as a consequence of the peak in the triplet part p d {ui)), G(0, V) 
decreases with applied bias voltage V for small V. However, when the window of width eV 
reaches the peak in the singlet contribution p s d (co), the average of the total spectral density 
p™{oj) is expected to increase and this leads to a peak in G(0, V) at finite bias voltages. 
This explains the result shown in Fig. [THat small temperatures. 

For other parameters, which would correspond to another experimental situation, in 
particular nearer to the transition, it might happen that p s d {uj) broadens as a consequence 
of the applied bias voltage and the structure with three peaks is absent. This is the case for 
the parameters of Fig. [12], where instead of three peaks, we obtain two small shoulders at 
positive and negative voltages, as shown by the full line in Fig. [121 The three-peak structure 
is however the general behavior expected for a system well inside the Kondo regime when 
E s — E t is larger that the width of the peaks in the spectral densities. 

Due to the difficulties in the calculation of nonequilibrium properties, an approximation 
usually made consists in taking the density of states calculated in equilibrium, and assume 
that it is constant with applied bias voltage. This approach, valid for non-interacting systems 
has been used for example in combination with first-principle calculations, to calculate the 



current through single molecules 68|]. Some results for the transport through two-level 



systems were obtained also in this approximation 



3l|, using a formula equivalent to Eq. 
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(|T7|) . but with the average density calculated for V = 0. 

From the result presented in Section IVlAI for finite V, it is clear that this procedure is not 
valid in general. An example is shown by the dashed line in Fig. O When the equilibrium 
densities are taken instead of the nonequilibrium ones, the effects of broadening of the 
spectral densities caused by the applied bias voltage are missed, and as a consequence, still 
the presence of two side peaks (although broad) is predicted for the conductance, in contrast 
to the full nonequilibrium calculation. Naturally, the splitting of the Kondo resonance at 
ligher bias voltages is also missed if frozen densities are assumed (see Section IVI Al and Ref 



M)- 



VII. SUMMARY AND DISCUSSION 

We have studied the singlet-triplet Anderson model (STAM) in which a configuration 
with a singlet and a triplet is mixed with another one with a doublet, via the hybridization 
with a conduction band. This is the simplest model that describes the conductance trough 
a multilevel quantum dot, hybridized with two leads in an effective one-channel fashion. As 
found earlier in intermediate valence systems p\ . the model has a quantum phase transition 
that separates a region with a singlet ground state from another one with a doublet ground 
state. This transition has been studied recently in transport measurements through Ceo 



quantum dots 391 ] . Our results provide an explanation to the observed behavior at both 
sides of the transition. In particular, in the region of parameters in which the ground state is 
a doublet, we obtain a zero bias conductance with a plateau at intermediate temperatures, 
which agrees with experiment. The three-peak structure observed in the nonequilibrium con- 
ductance as a function of applied bias voltage is also explained by the model. The separation 
of the electronic spectral density at the dot into two parts, which correspond to excitations 
involving either the singlet or the triplet, leads to a more transparent understanding of the 
underlying physics. In particular, the above mentioned plateau and the observed peak at 
finite bias voltages are due to singlet excitations. 

We have also studied several limits of the model, which allow us to shed light on the 
expected behavior of the conductance at very low temperatures and bias voltages (G(0, 0)), 
as the integer valence limits of the STAM, and an special quantum critical line in which the 
model can be mapped into an ordinary Anderson model plus a free spin 1/2. For V t = V s , 
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the latter model with additional exchange interaction H' [Eq. (|3T!) ] is equivalent to the 
STAM. An analysis of H' and the integer valent limits show that the value of G(0, 0) is 
consistent with a generalized Friedel-Luttinger sum rule derived recently 3l|, and which 
takes into account that the system is a singular Fermi liquid when the ground state is a 
doublet, in a similar way as the underscreened Kondo model j^. As it is clear form Eqs. 
[37]) and (1391) this means that an abrupt conductance change takes place at the transition 
3JJ . Although the NCA cannot reach zero temperature, our numerical results are consistent 



with this result. On the quantum critical line, the conductance is smoothly connected to 
that of the phase with a doublet ground state. 

Our NCA approach has the advantage over NRG that it can be rather easily extended to 
the nonequilibrium regime. When the configuration with an even number of particles is the 
favored one (as we have assumed here), it works very well on the quantum critical line and 
one expects it to be accurate enough on the "triplet" side of the transition (where the ground 
state is a doublet). However, on the "singlet" side, a spurious peak in the singlet part of the 
spectral density develops at very small temperatures and bias voltages, rendering our results 
quantitatively inaccurate for these parameters. We have found that similar difficulties arise 
when the configuration with an odd number of particles is favored, but now on the triplet 
side of the transition, including the quantum critical line. In fact, the ordinary Anderson 
model that comes out of the mapping described in Section fVAl has now an occupation near 
zero and develops a spurious peak at low temperatures and bias voltages. 

While our results on the triplet side of the transition agree with the experimental results 
of Roch et al. 39j, and the interpretation of them is rather simple, we cannot totally rule 
out the possibility that some of them are due to the NCA approximation, and that the 
physical explanation of the observed phenomena is different. For example, the plateau of 
the conductance on the triplet side of the transition has not been reported in previous NRG 



1 Q- 



calculations [3l|, |40(. We 
NRG at finite frequencies [' 
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ieve that this might be due to the lack of resolution of the 



64j This might be improved averaging over different shifted 



logarithmic discretizations (z averaging) 69|, [70j| and using recent developments (the full 
density matrix NRG) 7l[ Q] . In any case, taking into account the difficulties to extend the 



technique out of equilibrium [73j, it seems that a combination of both techniques (numerical 
renormalization group and non-crossing approximation) might be suitable to obtain further 
progress. 
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